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I. INTRODUCTION 


The nonlinear Schrodinger (NLS) equation plays an 
important role in many fields of nonlinear physics [l|-0. 
The cubic NLS equation is shown to be completely inte- 
grable for both self-focusing and defocusing Kerr nonlin- 
earities [bl-Q, which can be used to describe the propa¬ 
gation of optical pulses in Kerr-type optical media [3|-l5| , 
or to describe the dynamics of matter waves in Bose- 
Einstein condensates, known as the Gross-Pitaevskii 
equation [ikl. 

The NLS equations with real external potentials and 
gain-and-loss distributions have been studied in many 
works [mil because the refractive index of optical ma¬ 
terials is usually complex, i.e., n{x) = ur^x) + ini{x) 
with nR{x) and ni{x) being the real and imaginary parts, 
respectively. In optics, the propagation of a signal is 
stable unless the propagation constant of the light is 
in real spectrum range. This requirement can be ef¬ 
ficiently achieved if the gain-and-loss distributions in 
the medium are exactly balanced to ensure the relation 
n{x) = n*{—x) (or nR{x) = nR{—x) and ni{—x) = 
—ni{x)), which is known as the parity-time {VT) sym¬ 
metric systems [l8|. T^T-symmetry may exhibit en- 
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tirely real spectrum of the respective optical potential 
in some parameter regions, referred to as the unbroken 
PT-symmetry 18-2^. Beyond these regions, in the bro¬ 
ken T^T-symmetry, the spectrum becomes complex and 
propagating waves may be either grow or decay. 

Recently, various PT-symmetric potentials have been 
introduced to the NLS equations, which have been 
shown to possess stable and unstable solitons of differ¬ 
ent types [ill-lli, [sS-Sil. Examples include the NLS 
equations with T^T-symmetric Scarff II potential [lit 
251, periodic potential jl^, 2^ 271. h armonic poten- 
tial 1281, 1291, Rosen-Morse potential [SOl, Gaussian po- 


l3ll . l32| . sextic anharmonic double-well po- 
1, time-dependent harmonic-Gaussian poten- 
self-induced potential [35[, and etc. (see. 


tential 
tential 
tial [34 

e.g., |36l-l42| ). One the other hand, due to the signif¬ 
icant progress achieved in recent years on developing 
optical materials with adjustable refractive index, VT- 
symmetric optical systems made of solid-state waveg¬ 
uides and fiber networks j43l-[45j|. multi-level atomic sys¬ 
tems jibl-lioj. and microcavities iillii have been sug¬ 
gested or realized experimentally. These practical sys¬ 
tems provide a solid ground for the study of solitons in 
the NLS model with T^T-symmetric potentials. 

The T^T-symmetric potentials play a key role in the 
wave propagation in the NLS equations in fibre and 
waveguide optics. As an example, the wave propaga- 
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tion with a T^T-symmetric Scarff II potential in both 
linear and nonlinear models has been studied in recent 
years 


22l-l25|. The PT-symmetric Scarff II potential 


can support stable bright solitons in the NLS models 
within particular parameter regions More impor¬ 

tantly, the real and imaginary parts of PT-symmetric k- 
wavenumber Scarff II potential 5^ both approach to 
zero as \x\ oo [cf. Eq. (j4]), for Vq = 2, Wq = 1, k = \/2, 


we have V{x) ^ 3.25 x 10“^^ and W{x) ^ 1.04 x 10“^^ 
as X = 20], that is to say, they have limited effect on the 
nonlinear waves for the NLS equation only in a boundary 
region (e.g., \x\ < 20 for k = a/ 2) [cf. Eq. (|3j)] and can 
more easily support the existence of bright solitons in the 
NLS equation, compared with the T^T-symmetric har¬ 
monic and optical lattice potentials, which always have 
the effect on the nonlinear waves in the whole region (see, 
e.g, (22I. [2^. ) . However, a complete analysis of the 

soliton stability in a T^T-symmetric /c-wavenumber Scarff 
II potential is still lacking in a full parameter space (e.g., 
the space consists of amplitudes of real and imaginary 
parts of PT-symmetric /c-wavenumber Scarff II potential 
and the wavenumber /c, see Eq. (m). In addition, the 
study of optical solitons could be very useful in the infor¬ 
mation science. Eor example, one can use optical solitons 
in coding for secured optical communication, etc. 

In this work we propose a unified study of the opti¬ 
cal bright solitons governed by self-focusing and defo- 
cusing nonlinear Schrodinger (NLS) equations with gen¬ 
eralized PT-symmetric Scarff II potentials. Particu¬ 
larly, we consider a T^T-symmetric /c-wavenumber Scarff 
II potential and a multi-well Scarff II potential, re¬ 
spectively. Eor the /c-wavenumber Scarff II potential, 
we show that by using two /c-wavenumber rays and a 
/c-wavenumber parabola one can divide the parameter 
space {(Vb, ITo)|Vb > 0, Wq G M} (Vq and Wq are, respec¬ 
tively, the real and imaginary parts of the potential) in 
different regions, corresponding to unbroken and broken 
T^T-symmetry and the bright solitons with self-focusing 
and defocusing Kerr nonlinearities. Eor the multi-well 
Scarff II potential we obtain bright solitons by using a 
periodically space-modulated Kerr nonlinearity. Then, 
we analyze the linear stability of bright solitons for both 
potentials in details by using numerical simulations. Sur¬ 
prisingly, stable and unstable bright solitons are found 
in both regions of unbroken and broken PT-symmetry 
due to the existence of the nonlinearity. Eurt her more, 
we explore the bright solitons in three-dimensional (3D) 
self-focusing and defocusing NLS equations with a gen¬ 
eralized PT-symmetric Scarff II potential. 

Turning to the possible experimental implementation 
of the models we will study [see Eq. o with VT- 
symmetric /c-wavenumber and multi-well Scarff II poten¬ 
tials and m ], we can exploit a T^T-symmetric re¬ 
fractive index profile (e.g., the considered T^T-symmetric 
/c-wavenumber Scarff II potential) imprinted in a cold 
gas of two atomic isotopes of A-type configuration (say 


of ^^Rb and ^^Rb isotopes) loaded in an atomic cell, as 
suggested in Refs. [46145Q| . By the interference of two Ra¬ 
man resonances, the required spatial distribution of the 
refractive index (i.e., the /c-wavenumber and multi-well 
Scarff II potentials) can be achieved by a proper combina¬ 
tion of a control laser field and a far-off-resonance Stark 
laser field. Since the refractive index of the atomic vapor 
is determined by two external laser fields, whose intensi¬ 
ties can be lowered to several microwatts achievable for 
the lasers nowadays, the system has the advantages of 
real-time all-optical tunable capability and controllable 
T^T-symmetry accuracy. In addition, the studied model 
supports a large Kerr nonlinearity (which is at least 10^^- 
order larger than those measured for usual nonlinear op¬ 
tical materials) due to the Raman resonant character, 
which favors to the formation of optical bright solitons. 

The rest of this paper is arranged as follows. Sec. II 
gives the bright spatial solitons in one-dimensional 
T^T-symmetric potentials. Both T^T-symmetric k- 
wavenumber and multi-well Scarff II potentials are con¬ 
sidered. The linear stability of optical bright solitons is 
analyzed. Stable solitons are found for both unbroken 
and broken T^T-symmetric potentials. Sec. Ill gives the 
bright solitons in 3D self-focusing and defocusing NLS 
equations with a generalized T^T-symmetric Scarff II po¬ 
tential. Einally, the last section contains a summary of 
main results obtained in this work. 


II. BRIGHT SOLITONS IN 
ONE-DIMENSIONAL RT-SYMMETRIC 
POTENTIALS 

A. Nonlinear physical model and stationary 
solutions 

Now we uniformly investigate the exact localized solu¬ 
tion modes and their stabilities of self-focusing and de¬ 
focusing NLS equations with a PT-symmetric potential. 
In the dimensionless form, the physical model can be 
written as [ 2 ^ 

idt^ + dl'ijj + [V{x) -h iW{x)]'ip -h g\'ipf'ip = 0, (1) 

where dt = d/dt^ dx = djdx^ = '0(x,t) is a complex 
field of X, t {x and t are, respectively, dimensionless space 
and time), V{x) is a real external potential, W{x) is a 
real gain-and-loss distribution, g characterizes the self- 
focusing {g > 0) or defocusing {g < 0) Kerr nonlinearity, 
respectively. Eq. (HD is associated with a variational prin¬ 
ciple ^£('0)/^T* = 0 with the Lagrangian density 

dtp) = i {dd - 

+2[V{x)+iW{x)]\'ip\'^ + g\'ip\‘^. (2) 

The T^T-symmetric potential V{x) + iW{x) leads to the 
sufficient (but not necessary) condition V{x) = V{—x) 
and W{x) = —W{—x). The quasipower and power of 
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Eq. m are given by Q{t) = 2p{x,t)2p'^{—x,t)dx and 

P{t) = |'0(x, respectively. One can read¬ 

ily get that Qt = i gi;{x,t)^p*{-x,t)[\^{x,t)\‘^ - 

\^{—x,t)\‘^]dx and Pt = —2 j^'^W{x)\^jJ{x,t)\^dx. No¬ 
tice that if we set t ^ z in Eq. o, where 2; denotes 
the propagation distance [HI, then the following results 
about Eq. (HD still hold as t ^ z. 

We focus on the stationary solutions of Eq. ([T]) in the 
form '0(x,t) = where // is the real propaga¬ 

tion constant and the complex function (j){x) satisfies the 
stationary NLS equation with varying parameter modu¬ 
lation 

I^Hx) = + [V{x)+iW{x)](p{x)+g\(l){x)\'^(p{x), (3) 

which can be solved by using numerical methods. 

In the following, we consider two different types 
of T^T-symmetric potential, i.e., the T^T-symmetric k- 
wavenumber Scarff II potential and the periodically 
space-modulated Scarff II potential, and study the ex¬ 
act localized solutions and their linear stability in both 
regions of unbroken and broken T^T-symmetry. 


B. T^T-symmetric A;- wavenumber Scarff II potential 


To study the soliton solutions of Eq. m, we first con¬ 
sider the PT-symmetric /c-wavenumber Scarff II poten¬ 
tial [53I, , which is given as 

V{x) = Vo sech\kx ), W{x) = Wo sech(kx) tanh(/cx), (4) 


where k > 0 denotes the wavenumber, and Vb > 0 and 
Wo are real parameters and modulate amplitudes of the 
external potential and gain-and-loss distribution, respec¬ 
tively. The wavenumber k and amplitude Vb can mod¬ 
ulated the well width and depth of the potential E(x), 
respectively. Wo can modulate the impact of the gain- 
and-loss distribution. Eor the case /c = 1, T^T-symmetric 
A:-wavenumber Scarff II potential ffl becomes the usual 
T^T-symmetric Scarff II potential [l9|. Eigure[T] displays 
the T^T-symmetric A:-wavenumber Scarff II potential (|4)) 
for different amplitudes Vb, Wo and wavenumber k. 

It follows from Eq. (|4|) that E(x), W{x) ^0 as \x\ 

00 and/or k 00, and V{x) ^ Vb, W{x) 0 as x = 0 
and/or A: ^ 0. If VbWo 7^ 0, then V{x) and W{x) satisfy 
the family of elliptic curves: 


Vjx) 

Vo 


W^{x) 


1 

4 


(5) 


with the centers being (Vb/2,0) in the {V{x),W{x))- 
space . It is easy to see that the family of elliptic curves 
is dependent on the wavenumber k. 

In general, for the T^T-symmetric A:-wavenumber Scarff 
II potential S'pr = V{x) iW{x) with V(x) and W(x) 
given by Eq. 0 for any real parameters Vb, Wo there are 
the following three cases: 



FIG. 1. (color online) Real (a) and imaginary (b) parts of VT- 
symmetric A:-wavenumber Scarff II potential ( 0 . (a): Vo = 
3, A: = 1 /v^ (dashed-dotted line), Vo — 2^ k — \f2 (solid 
line), Vo = A; = 1 (dashed line), and (b): Wo — Z^k— l/\/2 
(dashed-dotted line). Wo — 2^ k — y/2 (solid line). Wo — k — 
1 (dashed line). 




FIG. 2 . (color online) The parameter space {(Vo,Wo)|Vo > 
0, Wo G M} can be divided in different regions by two families 
of A;-wavenumber open rays /i ,2 : Wo = ±(Vo -f A;^/4) with 
Vo > 0 and a family of A:-wavenumber parabola h : Wq = 
9A:^(Vo — 2 A:^), corresponding to unbroken and broken VT- 
symmetry and bright solitons with self-focusing (^ > 0 ) and 
defocusing ( 5 ^ < 0) Kerr nonlinearities. The tangent points of 
three curves are A ,2 = (17A;^/4, ±9A:^/2). 


i) When Wq = 0, the potential S'pr = V{x) is a real, 
Hermitian, and T^T-symmetric potential; 


ii) When Vb = 0, the potential S'pr = iW(x) is a pure 
imaginary, non-Hermitian, and PT-symmetric po¬ 
tential; 


hi) When VqWo 7^ 0, the potential S'pr = V{x) P 
iW{x) is a complex, non-Hermitian, and VT- 
symmetric potential. 


In the following, we mainly consider the potential 
S^pr = V{x) P iW{x) in Case hi) with Vq > 0. 
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TABLE L The regions for unbroken and broken PT-symmetry and the existence of bright solitons with self-focusing (p > 0) 
and defocusing {g < 0) nonlinearities in the parameter space {(Vb,Wo)|Vo > 0, Wo G M} [The signs ‘Yes’ and ‘No’ denote 
the corresponding problems exist and do not exist in the regions, respectively. The signs ‘+’ and ‘-’ stand for the union and 
difference of two sets, respectively]. 


Region 

Linear problem (^ = 0): 
unbroken PT-symmetry 

Linear problem (^ = 0): Nonlinear problem (^ > 0): 
broken PT-symmetry soliton (I10|) 

Nonlinear problem (^ < 0): 

soliton (ITOl) 

Ii + I2 + I3 + li 
Th — {^1,2} 

Yes 

No 

Yes 

No 

13 

Yes 

No 

No 

No 

II 

No 

Yes 

Yes 

No 

III 

Yes 

No 

No 

Yes 




FIG. 3. (color online) Real (solid lines) and imaginary (dashed lines) parts of bright soliton (ITOl) for (a) the self-focusing 
nonlinearity ^ = 1 with Vo = 5, Wo = 5.2, k — a/2 , and (b) defocusing nonlinearity g — —1 with Vo = 3, Wo = 0.5, k — 

(c) The transverse power-ffow vector (Poynting vector) S{x) given by Eq. (1161) for the self-focusing nonlinearity g = 1 with 
Vb = 5, Wo = 5.2, k — \f2 (solid line) and defocusing nonlinearity g = —1 with Vb = 3, Wo = 0.5, k = l/\/2 (dashed line). 


1. linear eigenvalue problem with unhrokon/hrokon 
VT -symmetry 

The linear eigenvalue problem 

L-ptHx) = Lvt = -dl - [V{x) + iW{x)] (6) 

with V(x) and W(x) given by Eq. (jUl exhibits an entirely 
real spectrum provided that two amplitudes Vb > 0, Wo 
and the wavenumber k satisfy 

\Wo\<Vo+’^, ( 7 ) 

where A and ff>(x) are eigenvalue and engenfunction, re¬ 
spectively. The inequality ([7)) is called the VT symmetry¬ 
unbreaking condition. If | Wo| > Vb + /c^/4, the spectrum 
becomes complex. For /c = 1, the conditio n 0 reduces to 
the well-known one, |Wo| < Vb + 1/4 (lol. [2^^ Therefore, 
two families of /c-wavenumber open rays li : 

Wo = Vo+k^lr ro>0 (8) 

and I 2 : 

(9) 


divide the parameter space {(Vb,Wo)|Vb > 0, Wo G M} 
into two regions, i.e., the unbroken T^T-symmetric region 
in/)I < Vb + A^^/4 and the broken PT-symmetric region 
iWoi >Vo^k^/4 (see Fig.EJ. 


2. The eonditions for the existenee of nonlinear modes 


For the given T^T-symmetric /c-wavenumber Scarff II 
potential (j4]) with Vb > 0 and Wo G M, Eq. admits 
the unified bright soliton for both self-focusing (^ > 0) 
and defocusing {g < 0) nonlinearities 

(fix) = W i - Vb + 2/^2^ sech(/cx)e*^^^\ (10) 


where= ^[Wq/( 9/c^) — Vb+ 2/c^] > 0, and the phase 
is 


Wo 

(fix) = —^ arctan[sinh(/cx)]. 

Orv 


( 11 ) 


In particular, when g = k = 1 or —g = k = the 
solution (ITOl) reduces to the known ones as given in 

Refs. SM. 


11^0 = -(Vo + fcV4), Vo > 0 
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Now we analyze the conditions for the existence of 
bright soliton m for parameters Wq, Vo, k, and g. For 
the positive (self-focusing) nonlinearity ^ > 0, the condi¬ 
tion for the existence of bright soliton ([TQ|) is given by 

>9k‘^{Vo-2k^). (12) 

The real and imaginary parts of the soliton are displayed 
in Fig. [3fa). For the negative (defocusing) nonlinearity 
^ < 0, the condition for the existence of bright soliton 
m is given by 

|Wo| < 3fcVFo - 2P, Vo > 2k^. (13) 

The real and imaginary parts of the soliton are illustrated 
in Fig. [3fb). 

Taking into account the conditions ([12]) and ([T3|) for 
the existence of bright soliton m for ^ > 0 and g < 
0, the parameter space {(Vb,lFo)|Vb > 0, IFo G M} can 
be divided into another two fundamental regions by the 
family of parabolas Is : Wq = 9k‘^{Vo — 2k‘^) with vertex 
being {2k‘^,0). These two conditions (fT^ and ([T3|) for 
bright solitons with ^ > 0 and ^ < 0, respectively, are 
shown in Fig. |2j 

Notice that the bright soliton with ^ > 0 (at¬ 
tractive case for Eq. O) can also exist for the region 
{(Vb,lFo)|Vb < 0, IFo G M} even if we do not study the 
7^T-unbroken/broken of the linear eigenvalue problem in 
this region. 

Two branches of /c-wavenumber open rays li ^2 given by 
Eqs. ([8|) snd (|9|) are tangent to the parabola Is : 

= 9k\Vo - 2k^) (14) 

at the points Pi ^2 = (17/c^/4, ±9/c^/2), respectively. In 
another word, |Vb + ^^/4| > 3/cVVb — 2k‘^ always holds 
except for two tangent points Pi, 2 , at which they are 
equal (see Eig. |2]and note [H^ for the details). 

Therefore, the parameter space {(Vb,lFo)|Vb > 
0, Wo G M} can be divided into eight regions, i.e., 

{(K),lFo)|Eo >0, VFo gM} 

(15) 

= II -M2 + I3 + II + III + /l + /2 + /S, 

where denotes the union of two sets, these regions 
Ii, 2,3 5 II? Ill denote open sets without the corresponding 
boundaries (see Eig. |2|). To make it become clearer, we 
summarize these different regions in Table [H One can 
see that the region for bright solitons with self-focusing 
nonlinearity contains both regions for unbroken and bro¬ 
ken PT-symmetry, whereas the region for bright solitons 
with defocusing nonlinearity contains only the region for 
unbroken PT-symmetry. 

Eor any complex solution 0(x) = \(l){x)\e^'^V) with ip{x) 
bering a real phase of 0(x), we find that its corresponding 
transverse power-flow or Poynting vector and the gradi¬ 
ent of phase obey the relation S{x) = f (00* — 0*0^,) = 


\(l){x)\‘^(Px. Eor the bright soliton ([TQ|) . its corresponding 
transverse power-flow or ‘Poynting vector’ is given by 

^(a:) =- Fo + 2 fc 2 ) sech^(fca;), ( 16 ) 

which implies that sgn(5'(x)) = sgn(IFo) for any x and 
^ = ±1 (see Eig. life)) since k > 0 and g{Wo /(9/c^) — Vb + 
2k^) > 0 are required for the conditions of the existence 
of bright solitons. Therefore, the power always flows in 
one direction, i.e., from the gain toward the loss. In 
addition, we have the conserved power related to solution 

dinjis P(i) = ^(g-Vo + 2fc2). 


3. Stability of nonlinear modes 


Next we study the linear stability of the bright soli- 
tons ([TQ|) in the above-mentioned different regions. To 
this end, we considered a perturbed bright soliton solu¬ 
tion [HU 


^(x, t) = |0(x)-h 6 F{x)e'^^^-\-G* {x)e | ( 




(17) 

where 0(x)e*^^ is a stationary solution of Eq. ([1]), e 1, 
and F{x) and G{x) are the eigenfunctions of the lin¬ 
earized eigenvalue problem. Substituting Eq. (ED in 
Eq. ^ and linearizing with respect to e, we obtain the 
following linear eigenvalue problem 


/ L g(l>^{x) 
\-9r\x) -L* 


F{x) 

G{x) 



,(18) 


where L = V{x) F iW(x) 2g\ 0(x) P — /i. 

The stability of the perturbed soliton fj{x,t) is related 
to the imaginary parts Im((5) of all eigenvalues S. If 
|Im((5)| > 0 , then the solution 'ip{x,t) will grow exponen¬ 
tially with t (i.e., it is unstable), otherwise the solution 
f^{x,t) is stable. In the following, we study the stabil¬ 
ity of exact soliton m under an initial random noise 
perturbation up to 2 % of its amplitude for different pa¬ 
rameters 

Since the regions Ii and I 3 are symmetric about the 
Vb-axes, we only need to consider one of them (say Ii). 
Thus we leave four regions for studying the linear sta¬ 
bility of the bright solitons (IIQ|. The regions Ii, I 2 and 
II correspond to bright solitons for the self-focusing non¬ 
linearity g = 1, whereas the region III corresponds to 
bright solitons for the defocusing nonlinearity g = —1. 
To be specific, we study the stability of bright solitons in 
the regions for unbroken and broken TT-symmetry with 
k = 1 (the usual TT-symmetric Scarff II potential) and 
k ^ 1. Notice that we only consider IFo > 0 as W{x) is 
an odd function. 

If k = 1 and g = 1: (i) for Vb = O-l, Wo = 0-5 in the 
region II (broken TT-symmetry), we numerically obtain 
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FIG. 4. (color online) The linear stability eigenvalues (left 
column) and propagation of soliton intensity with 

g = k = 1 (right column), (a)-(b) Vb = 0.1, Wo = 0.5 (in 
the region II), (c)-(d) Vb = 1.5, Wo = 0.2 (in the region I 2 ), 
(e)-(f) Vb = 5, Wo = 5.2 (in the region Ii), and (g)-(h) Vb = 
-0.1, Wo = 0.1 (in the region {(W, Wo)|Vo < 0, Wo G R}). 


the spectrum of eigenvalues 6 (see Fig. lU^a)) and a sta¬ 
ble propagation of the soliton intensity |' 0 (x,t)p shown 
in Fig. m^b); (ii) For Vb = 1-5, Wo = 0.2 in the region 
I 2 (unbroken T^T-symmetry), we numerically obtain the 
spectrum of eigenvalues S (see Fig. lU^c)) and a stable 
propagation of the soliton intensity shown in 

Fig.m^d); (hi) For Vb = 5, Wo = 5.2 in the region Ii (un¬ 
broken T^T-symmetry), we numerically obtain the spec¬ 
trum of eigenvalues 6 (see Fig. lU^e)) and an unstable 
propagation of the soliton intensity |' 0 (x,t)p shown in 
Fig. Iljf); (iv) For Vb = —0.1, Wo = 0.1 in the region 



FIG. 5. (color online) (a) The linear stability eigenvalues; 
(b) The stable propagation of soliton intensity with 

^ — 1. The other parameters are Vb = 3, Wo = 0.5, k = 1 
(in the region III). 


{(Vb, Wo)|Vb < 0 , Wo G M}, we numerically obtain the 
spectrum of eigenvalues 6 (see Fig.|4{g)) and an unstable 
propagation of the soliton intensity |' 0 (x,t)p shown in 
Fig. |l{h). If k = 1 and g = —1, for Vb = 3, Wo = 0.5 
in the region III (unbroken T^T-symmetry), we numeri¬ 
cally obtain the spectrum of eigenvalues S (see Fig.[5{a)) 
and a stable propagation of the soliton intensity |' 0 (x, t) p 
shown in Fig. [5{b). We stress that the stability (insta¬ 
bility) of the bright solitons in broken (unbroken) PT- 
symmetry region is due to the existence of Kerr nonlin¬ 
earity. 

If k ^ 1 and g = 1: (i) for Vb = 0.05, Wo = 

0.25, k = 1/a/ 2 in the region II (broken T^T-symmetry), 
we numerically obtain the spectrum of eigenvalues 6 
(see Fig. [ 6 {a)) and a stable propagation of the soliton 
intensity |'0(x,t)p shown in Fig. [ 6 {b); (ii) For Vb = 
2.16, Wo = 0 . 12 , k = \/T^ in the region I 2 (unbroken 
T^T-symmetry), we numerically obtain the spectrum of 
eigenvalues S (see Fig. [ 6 {c)) and a stable propagation of 
the soliton intensity |'0(x,t)p shown in Fig. [ 6 {d); (hi) 
For Vb = 10 , Wo = 10.4, k = \/2 in the region Ii (un¬ 
broken PT-symmetry), we numerically obtain the spec¬ 
trum of eigenvalues h (see Fig. [ 6 {e)) and an unstable 
propagation of the soliton intensity |' 0 (x,t)p shown in 
Fig. [ 6 {f); (iv) For Vb = —0.05, Wo = 0.05, k = l/\/2, 
in the region {(Vb,Wo)|Vb < 0, Wo G M}, we numeri¬ 
cally obtain the spectrum of eigenvalues 6 (see Fig.[ 6 {g)) 
and an unstable propagation of the soliton intensity 
|'0(x,t)p shown in Fig. [ 6 {h). If /c 7 ^ 1 and g = — 1, 
for Vb = 6 , Wo = 1, k = \/2 in the region III (unbroken 
T^T-symmetry), we see that |Im(h)| = 0 (see Fig. [3^a)) 
and a stable propagation of the soliton intensity 
is obtained as shown in Fig. [7{b). 


C. 'PT-symmetric multi-well ScarfF II potential 


Now we consider another type of PT-symmetric Scarff 
II potential, i.e. the PT-symmetric multi-well Scarff II 
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FIG. 7. (color online) (a) The linear stability eigenvalues; 
(b) The stable propagation of soliton intensity |'0(x,t)|^ with 
g = —1. The other parameters are Vo = 6, Wo = 1, /c = a/2 
(in the region III). 



FIG. 6. (color online) The linear stability eigenvalues (left 
column) and propagation of soliton intensity |'0(x,t)|^ with 
g = I (right column), (a)-(b) Vo = 0.05, Wo = 0.25, k = 
1/a/ 2 (in the region II), (c)-(d), Vo = 2.16, Wo = 0 .12, k = 
(in the region I2), (e)-(f) W = 10 , Wo = 10.4, k = 

(in the region Ii), and (g)-(h) Vo = —0.05, Wo = 0.05, k = 
I/V 2 (in the region {{Vo, Wo)\Vo < 0, Wo G R}. 


potential, reading 


V{x) 




p + 2 — CT cos(co’x) sech^(x), 


(19) 


FIG. 8. (color online) The external potential —V(x) exhibit¬ 
ing the multi-well structure for (a) cr = 2 and (b) a = —2. 
Other parameters are Wo = 0.2, cu = 1,6. 


(a = —2, cj = 1), and multi-well (a = ±2, cj = 6) struc¬ 
tures, respectively (see Fig. [8]). In fact, for the fixed pa¬ 
rameters Wo and cu, one can change the parameter a to 
control the number of the wells. Notice that in the figure 
we plot the profile of —V{x) because the Hamiltonian of 
Eq. ^ is —dl — V{x) — iW{x). As cj ^ 0 (or cj ^ 1), 
the linear problem associated with the T^T-symmetric 
potential dm admits an entirely real spectra provided 
that I Wo I < VEo /9 + 9/4 — cr. Ifcj^l, the condition for 
the unbroken-symmetry is complicated. For the fixed 
parameter Wo = 0.2, we find the regions of the unbroken- 
VT symmetry: (i) uj = 0.5, —280 < cr < 2.5; (ii) cj = 1, 
—53.3 < (T < 3.17; (hi) uu = 2, —13.1 < cr < 6.65; and 
(iv) cj = 5, -13.9 < cr < 11.33. 

For the PT-symmetric multi-well Scarff II potential 
m, we consider the periodically space-modulated non¬ 
linearity g ^ g{x) in Eq. (|3j) with the form 

g{x) = crcos(co’x). (20) 



W{x) = Wosech(x) tanh(x), 

where Wo, a G M and cj > 0 denotes the wavenumber. If 
UJ = 0 and a = Wq / 9 + 2 — Vb, then V{x) becomes the 
usual Scarff II potential which exhibits the shape 

of single well, but for non-zero wavenumber uj, V{x) ex¬ 
hibits the shape of multi-well, which can provide more 
abundant structures. Eor Wo = 0.2, the external poten¬ 
tial V{x) exhibits double-well {a = 2, uj = 1), singe-well 


Eor the zero wavenumber uj = 0, nonlinearity g{x) is 
a constant, g{x) = cr, however, for the given non-zero 
wavenumber uj, the nonlinearity g{x) is either positive 
or negative as x increases. Based on the potential ([T9]) 
and the nonlinearity dsoj we have exact bright solitons 
of Eq. (|3|) 

0(x) = sech(x)e*^^^\ 


( 21 ) 
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where /i = 1 and the phase is 
Wo 

ip{x) = -arctan[sinh(x)]. (22) 

O 

The corresponding transverse power-flow or Poynting 
vector is given by 


S{x) = -(p*(px) = i^sech^(a;), (23) 

which implies that sgn{S{x)) = sgn(]To) for any x. 
Therefore, the power always flows in one direction, i.e., 
from the gain toward the loss. Moreover, the conserved 
power is P{t) = = 2. 

Now we consider the linear stability of the soliton so¬ 
lution (ED using Eq. m with g being replaced by the 
periodic function dsni) with Wo =0.2 and different val¬ 
ues of cj, a. We find that for the fixed oo the soliton m 
is more (less) stable if cr < 0 (cr > 0). For examples, we 
show the linear stability eigenvalues and the stable prop¬ 
agation of the soliton ([21]) for cj = 1, cr = 2 (see Fig. Ha)- 
(b)), cj = 1, cr = —2 (see Fig. |9]^c)-(d)), uj = 0.5, cr = —3 
(see Fig.|9]^e)-(f)), and cj = 0.5, cr = 1 (see Fig.|9]^g)-(h)). 

For the fixed parameters Wo = 0.2, cj = 0.5, when 
cr = —310, for which the T^T-symmetry is unbroken, the 
soliton is still stable (see Fig.flOfb)). however, when 
cr = 2.5, in which T^T-symmetry is unbroken, the soliton 
m becomes unstable (see Fig. flOf d)). Thus, for Wo = 
0.2, UJ = 0.5, the parameter cr has significant influence on 
the linear stability of the soliton solution (|2T]) . A negative 
cr favors to obtain the stable solitons. 


III. SOLITONS IN THREE-DIMENSIONAL 
PT-SYMMETRIC POTENTIAL 



0 


-20 -10 0 10 20 
Re (5) 




Re(6) 


Intensity 



-5 0 5 


X 


Intensity 



-5 0 5 


X 


We now consider the 3D NFS equation with the VT- 
symmetric potential 

idti^+\7^'ip+[V{x,y,z)+iW{x,y,z)]'ip+g\'ip\^'ip = 0,(24) 

where t/; = ^{x^y^z^t) is a complex field with respect 
to x,y,z,t G M, = 5^ + 9^ + 9^, V{x,y,z) and 
W{x^y^z) are both real-valued functions related to ex¬ 
ternal potential and gain-and-loss distribution, respec¬ 
tively, and is a real constant with g = ±1. The 
T^T-symmetric potential requires the sufficient (but not 
necessary) condition V{x^y^z) = V{—x^—y^—z) and 
W{x^y^z) = —W{—x^—y^—z). The stationary solu¬ 
tion can be solved in the form '0(x, y, z) = 0(x, y, z)e^^^^ 
where g is the real propagation constant and 0(x, y, z) 
satisfies 


FIG. 9. (color online) The linear stability eigenvalues (left 
column) and stable propagation of soliton intensity 
(right column), (a)-(b) Wo = 0.2, cj = 1, cr = 2; (c)-(d) 
IFo = 0.2, a; = 1, cr = -2; (e)-(f) Wo = 0.2, a; = 0.5, cr = -3; 
(g)-(h) Wo = 0.2, UJ = 0.5, cr = 1. All panels correspond to 
the unbroken PT-symmetry. 


We consider the generalized 3D T^T-symmetric Scarff 
II potential 


V{x,y,z) 


W{x,y,z) 


= F (^ + 2fc^)sech2(fc^7/) 

r]=x,y,z ^ V / 

= Wo J2ri=x,v,z sech{krjT]) tanh(fc^r?), 


(26) 


IJ(t) = V^<()-|- [V{x,y,z) + iW{x,y,z)](l) +(25) 


where kr^ > 0 {g = x^y^ z) are the wavenumbers m x^y^z- 
directions, respectively, and Wq, are real constants. 
In particular, the T^T-symmetric potential (|26|) reduces 
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FIG. 10. (color online) The linear stability eigenvalues (left 
column) and propagation of soliton intensity |'0(x,t)|^ ()21D 
(right column), (a)-(b) a — —310, corresponding to the bro¬ 
ken PT-symmetry; (c)-(d) <j — 2.5, corresponding to the un¬ 
broken PT-symmetry. Other parameters are Wq = 0.2 and 
a; = 0.5. 


to the Scarff II potential (j4]) when rj = x and = 1. 
Since x,y,z are symmetric in V{x,y,z) and W{x,y,z) 
thus we plot their profiles in (x, ^)-space with z = 0 (see 
Fig. [TTj). For the self-focusing nonlinearity > 0, the 
external potential V{x,y,0) exhibits the different profiles 
as g becomes large (see Fig. [TT]fa)-(c)). 

For the above-mentioned 3D potential ([26|) we obtain 
the exact bright solitons of Eq. (ESJ 

0(x, y, z) = (l)Qsech{kxx)sech{kyy)sech{kzz)e'^'^^^''^'^\{27) 

where y = k‘^ ^ ky + k^ and the phase is 

(p{x^y^z) = ^^ ^ arctan[sinh(/c^7^)]. (28) 

o 

ri=x,y,z 

The real and imaginary parts and intensity of the solution 
(j27|) are shown in Fig.[T^a)-(c) for 0o = 1 and Wo = 0.5. 

The velocity field v(x, y, z) of the solitons (|27|) have 
the form 

v = Vip{x,y,z) 

(29) 

= ^[k^^sech.{kxx), ky^sech.{kyy), k^^sech.{kzz)], 

which is shown in Fig. [T2ja) with Wo = 0.5. It fol¬ 
low from Eq. (|29)) that the divergence of velocity field 
v{x,y,z) (alias the flux density) is given by 

divv(a;, y, 2:) = V-v{x,y,z) = V'^(p(x,y,z) 

(30) 

= T,y=x,y,z sechikyv) tanh(^7?), 



FIG. 11. (color online) The external potential V{x,y,0) with 
(a) 9 = 1, (b) £/ = 3, (c) g = 8, and (d) g = -1; (e) The gain- 
and-loss distribution W(x,y,0). The parameters are Wo = 



FIG. 12. (color online) (a) the velocity field v{x,y,z) (1^ : 
(b) the flux density divv(x,y,0) (l3Ql) . The parameter are 
Wo = 0.5 and kx = ky = kz = 1. 


which measures the flux per unit area and is dependent 
on both Wo and space position. Fig.[T2fb) shows the flux 
density in (x, y, 0)-space. In addition, we also have the 
relation 

divv(a;, y, 2 ) = V^(f{x,y,z) = - ^W{x,y,z). (31) 

From Eq. (l3Q|) we have the following proposition. Eor 
the given parameter Wo > 0, we have: 

i) the fluid flows outward if f{x,y^z) < 0; 




















10 





FIG. 13. (color online) The (a) real part, (b) imaginary part, 
and (c) the intensity of the soliton solution (1271) : (d) The 
transverse power-flow vector (Poynting vector) S{x, y, z). The 
parameter are Wo = 0.5 and kx — ky — kz — 



ii) the fluid flows inward if f{x,y,z) > 0; 

iii) the fluid does not flow if f{x^y^z) = 0. 

where we have introduced f{x,y,z) = 

y ^sech{kr^y)ta,nh{kr^r]). For the case Wo < 0, we 
also have the corresponding results. 

For any complex solution (j){x^y^z) = 

|0(x, we find that its corresponding 

transverse power-flow or Poynting vector and the 
gradient of phase obey the relation 

S{x,y,z) 

(32) 

= \(t>{x,y,z)\‘^V(f{x,y,z). 

The transverse power-flow or Poynting vector related to 
the solution (j27j) is given by 


S{x,y,z) = 4>l sech^ v(a;, y, z) 


r]=x,y,z 
2 




sech (k^ri) 


(33) 


3 ^^r]=x,y,z 

X [k~^sech{kxx)^ ky^sech{kyy)^ k~^sech{kzz)] 

for either g = 1 or —1 (i.e. it does not depend on the sign 


of the nonlinearity), which is exhibited in Fig. fTdT d). In 
addition, the conserved power related to the solution 
(|27|) is given by P{t) = S(j)Q/{kxkykz) for either ^ = 1 or 

g = -!• 

IV. CONCLUSION AND DISCUSSION 


In conclusion, we have presented a unified theoret¬ 
ical study of the optical bright solitons governed by 
self-focusing and defocusing NLS equations with VT- 
symmetric Scarff-like II potentials. Particularly, a VT- 
symmetric /c-wavenumber Scarff II potential and a multi¬ 
well Scarff II potential are considered, respectively. For 
the /c-wavenumber Scarff II potential, the parameter 
space can be divided into different regions, correspond¬ 
ing to unbroken and broken T^T-symmetry and bright 
solitons with self-focusing and defocusing Kerr nonlin¬ 
earities. For the multi-well Scarff II potential the bright 
solitons can be obtained by using a periodically space- 
modulated Kerr nonlinearity. The linear stability of 
bright solitons with T^T-symmetric /c-wavenumber and 
multi-well Scarff II potentials is analyzed in details by 
using numerical simulations. Stable and unstable bright 
solitons are found in both regions of unbroken and bro¬ 
ken PT-symmetry due to the existence of the nonlinear¬ 
ity. Furthermore, the bright solitons in 3D self-focusing 
and defocusing NLS equations with a generalized VT- 
symmetric Scarff II potential are explored. The used 
method can also be used to study optcial solitons in other 
T^T-symmetric k-wavenumber potentials. 

The results we obtained in this work provide a new 
way of control over soliton stability by using generalized 
T^T-symmetric Scarff II potentials in different parameter 
domains. This may have potential applications in the 
field of optical information transmission and processing 
based on optical solitons in nonlinear dissipative but VT- 
symmetric systems. 
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